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Abstract 

The coefficient of of a formal Laurent series f(x) is called the formal residue of 
| f(x). Many combinatorial numbers can be represented by the formal residues of hyper- 

geometric terms. With these representations and the extended Zeilberger's algorithm, we 
generate recurrence relations for summations involving combinatorial sequences such as 
Stirling numbers and their g-analogue. As examples, we give computer proofs of several 
known identities and derive some new identities. The applicability of this method is also 
'— studied. 
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q '. 1 Introduction 

Finding recurrence relations for summations is a key step in computer proofs of combinatorial 
.^J | identities. In 1990's, Wilf and Zeilberger [20,21] developed the method of creative telescoping 

^ ■ to generate recurrence relations for hypergeometric summations. Since then, many extensions 

and new algorithms have been discovered and designed for various kinds of summations. See, 
for example, [1,2] for holonomic sequences, [3,19] for multivariable hypergeometric terms, [17] 
for nested sums and products, [11, 12] for Stirling-like numbers, [5, 13, 14] for non-holonomic 
sequences. 

Our approach is motivated by the work of Chen and Sun [4]. By using the Cauchy contour 
integral representations, they transformed sums involving Bernoulli numbers into hypergeo- 
metric summations. Then the recurrence relations for the sums can be derived by the extended 
Zeilberger's algorithm [3]. 

In the present paper, we combine the formal residue operator and the extended Zeilberger's 
algorithm to generate recurrence relations for combinatorial sums. With this residue method, 
we give computer proofs of some known identities and derive some new identities. Moreover, 
we study the applicability of this method. We show that in the case of one variable, it is 
equivalent to the Sister-Celine's method. 

We note that Egorychev [7] provided integral representations for many combinatorial num- 
bers and used them to prove combinatorial identities. Fiirst [8] reformulated Egorychev's 
method in terms of formal residue operators. Egorychev transformed sums into geometric 
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series and then evaluated them by some manipulation rules. We transform sums into hyperge- 
ometric sums and find recurrence relations they satisfied. 

The paper is organized as follows. In Section 2, we describe the residue method. Then 
in Section 3, we give several examples involving Stirling numbers of both kinds. Section 4 is 
devoted to deriving two new Stirling number identities. In section 5, we consider the g-Stirling 
numbers as well as other combinatorial sequences which also fall in the scope of our method. 
Finally in Section 6, we study the applicability of the residue method. 



2 The method of residue 

Let K be a field and K((z)) be the set of formal Laurent series in the indeterminate z over K. 
For any element 

oo 

f(z) = a nz n € K((z)), (2.1) 

n=riQ 

the formal residue operator res (or, res if no confusion) is defined by 

z 

res f(z) = res f(z) = a_i. 

z 

Clearly, the k-th coefficient of f(z) can be represented by formal residue as follows 

/(*) 

We see that this representation is equivalent to the Cauchy integral representation of a*,, 

1 I /(*) , 

Based on the formal residue, we give a computer assisted method to derive recurrence 
relations for sums involving non-hypergeometric sequences. Consider a definite sum with the 
form of 



f(n)= J2 F ( n ' fc )> 



k=— oo 

where n = (m, . . . , n r ) is the vector of parameters. The residue method consists of the following 
three steps. 

1. Rewrite the summand F(n,k) as res F(n,k, z), where F(n,k,z) is a hyper geometric 

z 

term. 

2. Take a finite subset 5cN r and apply the extended Zeilberger's algorithm to the similar 
terms {F(n + a, k, z)} ae s, where N denotes the set of non-negative integers. If the 
algorithm succeeds, we thus obtain a relation of the form 

PUn)F(n + a, k, z) = A k G(n, k, z), (2.2) 

aeS 

where p a {n) are polynomial coefficients independent of k and z and G(n,k,z) is a hy- 
pergeometric term similar to F(n,k,z). 
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3. Summing over k and applying the operator res^, we are led to a recurrence relation for 
the sum f(n), 

/ Pa(n)f(n + ot) = resG(n, +00, z) — resG(n, —00, z). 

— J z z 

Remark. In most cases, G(n, k, z) is finite supported and hence we do not need to calculate 
res z G(n, +00, z) and res 2 G(n, —00, z). 

To conclude this section, we give an example to illustrate the method of residue. More 
examples can be found in Sections 3-5. 

Example 1. We have [9, identity (6.15)] 

Y,(fys 2 (k,m) = S 2 (n + l,m + l), (2.3) 

where S 2 (n, m) is the Stirling number of the second kind. 
Proof. It is known that 

00 k 
£ 52(n ' = (l-z)(l-2z)...(l-kzY 



n=k 

Therefore, 



z k 



S ^ k) = T z^{l-z){l-2z)...{l-kzY (2 ' 4) 



Denote the left hand side of (2.3) by L(n,m). We thus have 

L(n, m) = res ( " ) , , , — . 



k 

Now consider the inner summand 

C(n, m, k, z) ■ 



M ^ +i n™i(i-^)' 

Applying the extended Zeilberger's algorithm to the four similar terms 

C(n + i,m + j,k, z), < i,j < 1, 

we find that 

C(n + 1, m + 1, k, z) — (m + 2)C(n, m + 1, k, z) — C(n, m, k, z) = A^- 



(n + 1 - fc)(l - (m + l)z)' 
Summing over k and applying the formal residue operator, we derive that 

— L(n, m) — (m + 2)L(n, m + 1) + £(n + 1, m + 1) = 0. 

This agrees with the recurrence relation satisfied by ^(ra + l,m + 1)- Finally, the identity 
follows by checking the initial values 

L(0,m) = S 2 (0,m) = S 2 (l,m + l), L(n, 0) = 1 = S 2 (n + 1, 1). I 

We remark that most of the sums appearing in this paper can also be treated by Koutschan's 
implementation of the creative telescoping algorithm on non-holonomic sequence (for more 
detail, see [13]). The only exception is Example 2. 
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3 Stirling number identities 



In this section we shall provide several examples involving Stirling numbers of both kinds to 
illustrate the residue method. Recall that 

n 

5i(n, k)z k = (z)n = z(z-l)---(z-n + l), 

k=0 

and 



J2S2(n,k)z n 



z k 



(1- z)(l - 2z) ■ ■ ■ (1 - kzV 

n=k 

where S\(n,k) and 52 (n, k) are Stirling numbers of the first kind and of the second kind, 
respectively. We thus have 

5i(n, k) = res ■ 

and 



z Z^~^~^ 



S2 (n, k) = res 



z k 



z z n+1 {\ - z){\ - 2z) ■ ■ • (1 - kz\ 

It is worth mentioning that we use the ordinary generating functions of Stirling numbers 
instead of their exponential generating functions, which have been extensively used in [8]. Let 



k 



Fi (n, k, z) = -Trpr and F 2 (n,k,z 



J+ 1 zy ' ' ' z n + 1 (l-z)( y l-2z)---( y l-kz)' 

We see that both Fi(n, k, z) and F2(n, k, z) are hyper geometric terms of n and k. Let TV" and K 
be the shift operators with respect to n and k, respectively. Denote the ring of linear difference 
operators with rational coefficients by 

K(n,k)(N,K) = J ^^njfatylSrK' : I, J G N, r i:j (n,k) G K(n, k) 1 . 

[ i=0 j=0 J 

We see also that 

Ann5i(n, k) C AnnFi(n, k, z) and Ann52(n, k) C AnnF2(n, k, z), 

where 

Ann/(n,A;) = {L G K(n, fc)(iV, : Lf(n,k) = 0}. 

Given a function F(n,k), we denote by F(n,k,z) the function obtained from F(n,k) by 
replacing 5i(n, fe) and S2(n,k) with Fi(n,k,z) and ^(n, /c,z), respectively. Suppose that 
there exist Q G K(n, /c)(N, and L G K(n)(N) such that 

L - (K - l)Q G AnnF(n,/e). 

Then we also have 

L - (K - l)Q G Ann F(n, /c, z), 

which leads to an equation of the form (2.2). The extended Zeilberger's algorithm will succeed 
in finding such L and Q. This fact indicates that the residue method always works as long as 
the existence of such L and Q is guaranteed. 

With the residue method, we can prove all identities on Stirling numbers appeared in [11]. 
Moreover, we can deal with sums involving products of Stirling numbers, typically are identities 
(6.24), (6.25), (6.28) and (6.29) in [9]. Here we only give two examples. 
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Example 2. We have [9, identity (6.24)] 

Y,Si(k,m)S 2 (n + l,k + l) = ( H \ (3.1) 

k 

Proof. Denote the left hand side by L(n,m). We have 

x k+1 (y)k 



L(n, m) = res res > 

X 11 J 



y ^ x n + 2 \[ k +l{l -ix)y m+1 ' 
For the inner summand F(n,m, k), Gosper's algorithm gives 

F(n, m, k) = G(n, m, k + 1) — G(n, m, k), 

where 

x k i.y^)k i 

G( n. in. k) — 



- n+1 nti(l - ix) V m+1 l-x{l + y)' 
Since the denominator contains l — x(l + y) as a factor, we are unable to deduce a closed form 
of res x iesy G. However, summing over k from to n, we get 

T , x if 1 (i/)n+l 1 \ 

L(n,m) = res res-- — — -— tt- — -== + 



^ y (l-x(l + y))\ Yli^ii 1 -^) y m a^+i y m+i^ 

Notice that 

1 1 (v)n+i . _ 

res res ; r- n = = res (J = (J, 

y x (l- x (l + y))Yl^(l- ix ) y™ y 

and 

1 1 (l + y) n (n 

y x (1 - x(l + y)) x" +1 y m+1 2/ y m+1 \m 
This completes the proof. I 

Example 3. We have [9, identity (6.28)] 

£ Q 5 2 (fc, 05 2 (n - k,m) = ^ + m ) 5 2 (n, / + m). (3.2) 

Proof. Denote the left hand side by L(n,m,l). We have 

t / n ( n \ xl y m 

Lin, m, I) = res res > ; . 

- V Y \kj x k+l y n-k + l n j =i(1 _ ix) n m _ jy) 

For the inner summand F(n,m,l, k), the extended Zeilberger's algorithm gives 

-F(n, m+1, l)-F(n, m, l+l)-(m+2+l)F(n, m+1, /+1)+F(n+1, m+1, Z+l) = A k G(n, m, I, k), 

where 

/V\ -kx l y m 

G{n, m,l, k) - 



kj( n + l- A;)x*+ V" fc+1 U l +\(1 - ix ) Yl™+\1 - jy) ' 
Summing over k and applying the operators res^ and res,,, we get a recurrence relation 

-L(n, m + 1, /) - L(n, m, I + 1) - (m + 2 + Z)L(ra, m + 1, / + 1) + L(n + 1, m + 1, / + 1) = 0. 

It is easy to check that the right hand side of (3.2) satisfies the same recurrence relation. 
Finally, the identities holds by checking the initial values 

L(0, m, /) = <5 ,/<5o,m, L ( n i 0, = S 2 (n, I), L(n, m, 0) = S 2 (n, m). 
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4 New identities 



In this section, we use two examples to illustrate how to discover new identities by the residue 
method. In the first example, we generate new identities by introducing a new parameter in the 
original summand. While in the second example, we use Zeilberger's algorithm to construct 
new identities, as done by Chen and Sun [4]. 

We first consider the identity 

k=0 

which was proposed by Kauers and Sheng-Lang Ko as the American Mathematical Monthly 
Problem 11545. It was proved by Fiirst [8] by the residue representation 

n\ ln fc (l + z) 
Si(n,k) = — res- 



kl z z n+1 
In fact, this identity can be generalized as follows. 

Theorem 4. Let n and m be nonnegative integers. Then we have 

A , 2n \ n/ , f(2n-l)!!, if m = 0, 

> -If , ]S 1 (n + k,k + m ) = < v 4.2 
J-^ \n + kj K I 0, ifm>\. 

k=—m 1 ' j — 



Proof. Denote the left hand side of (4.2) by L(n,m). By the residue representation, we have 

2n \ (z) n +k 



L ( „,„)= r? E(-i)'(„ 2 ; fc ) 



k=—m 

The extended Zeilberger's algorithm gives the recurrence relation 

2(n + l)(2n + 3)L(n, m) - ( 2n + 3 )( 4n + 3 ) L(n + i, m ) + L {n + 2, m) 

+ 2(n + l)(2n + 3)L(n + l,m + l) = 0. (4.3) 

We now prove that L(n, n — r + 1) = for n > r by induction on the non-negative integer 
r. Since S\(n + k,n + 1 + k) = for any integer k, we have L(n, n + 1) = 0, i.e., the assertion 
holds for r = 0. For r = 1, we have 



L(n,n)= E(-l) fc ( n 2 + \)=0, n>l. 



k=—n 

Now suppose that the assertion holds for 1 < r < ro where ro > 1. The recurrence relation 
(4.3) implies that 

(2n — l)(An — 5) 

2(n - l)(2n - l)L(n - 2, n - r ) - n ^L(n - 1, n - r ) + L(n, n - r ) 

+ 2(n - l)(2n - l)L(n - 1, n - r + 1) = 0. 

By induction, we have 

L(n — 2, n — ro) = L(n — 1, n — ro) = L(n — 1, n — ro + 1) = 0. 
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Therefore, L(n,n — ro) = 0, which completes the induction. Notice that the assertion is 
equivalent to the statement L(n,m) = for any non- negative integers n and m > 1. 

For m = 0, the recurrence relation (4.3) becomes 

2(n + l)(2n + 3)L(n, 0) - ( 2n + 3 )( 4 ^ + 3 ) L(n + 1,0) + L{n + 2, 0) = 0. 

It is easy to check that (2n — 1)!! satisfies this recurrence relation and coincides with the initial 
values L(0,0) = L(1,0) = 1. I 

R. Sitgreaves [18] found the following identity (see also [7]). 



n / i \ 

^2 ( , J(-l) k S 2 (n + m - k,n - k) = 0, m > 0, n>m + l. 

z — n V / 



(4.4) 



fc=0 

From this result, we can establish the following theorem. 



Theorem 5. For nonnegative integers n > m > 0, we have 

yv / n + m + 1\ ( _ 1)fe)S2(n + m _ fc) n _ fc j = (_i)n+m m ,_ (4 5) 

Proof. Denoting the left hand side of (4.4) by L(n,m), we have 

k 



Lin m) = res^( n + m )(-l) k - 

z ^ V *: / z n+m-fe + i n n-fe (1 



For the inner summand F(n, m, /c), the original Zeilberger's algorithm gives 

(n + m + l)F(n, m, fc) - z(m + l)F(n, m + 1, fc) — zF(n, m + 2, fc) 

(n + m + l)kF(n, m, k) 



A* 



(n + m + 1 — /c)(n + m + 2 — /c)z 
Summing over fc and applying the residue operator, we obtain 

^ J (— l) fe S , 2(n + m — k, n — k) 

I J(-l) fe 5 2 (n + m + l-A:,n-A:) = 0. 

fe=o ^ ' 

Denote the left hand side of (4.5) by S(n,m). Substituting L(n, m) = in the above identity, 
we deduce that 

(m + l)5(n, m) + S(n, m + 1) = 0. 

Thus we have 

S(n, m + 1) = (-l) m+1 (m + l)!5(n, 0). 

Note that 

^ °) = e ( n r) ( - i)fc ^ (n - n - fe) = e ( n r) = ( - i)n - 

we finally derive that 

5(n,m + l) = (-l) n+m+1 (m + l)!, 
as desired. I 



5 More combinatorial sequences 



It is readily seen that our approach is also appliable to many other combinatorial sequences as 
long as the corresponding generating function is hypergeometric. More generally, the residue 
operator can be replaced by any linear operator L. For example, a classical treatment for 
identities involving harmonic numbers H n = Ylk=i ( see [15]) is to use the fact 



H n = 5 



n + x 

x 



where Sf(x) = ^\ x=0 . 

Here we list several sequences which could be treated by this method. 
g-Stirling numbers 

A kind of g-analogue of Stirling numbers is given by [6, 10] 

Sf( n , k) = S q (n - 1, k - 1) - [n - l]S*(n - 1, k), Sf(0, k) = 5 , k , 

S q (n, k) = S q 2 {n - 1, k - 1) + [k]S q 2 {n -l,k), S q (0, k) = d , k , 

where 

V-jT =1+ + ... + n-l_ 
1-q 

Their generating functions are 

n n—1 

£ Sl(n, k)z k = J] (z - [k]), £ S q (n, k)z n = . 

k=0 k=0 n>k fl (1 - [l]z) 

i=l 

Thus we have 

71-1 



n(*-H) 

i=0 QI l \ Z 



k—n—l 



S 9 (n, k) = res t— t , S 2 (n, k) = res 



Note that we also have 



k+1 , — k 

n (i-[<» 
i=i 



n-l 



Ann5f(n, fc) C Ann - and Ann S 2 (n, k) C Ann . 



n 

1=1 



Using these representations and the g-analogue of the extended Zeilberger's algorithm, we 
can derive recurrence relations for sums involving g-Stirling numbers. For instance, let us 
consider the sum (see [12]) 

n 

L(n,m) = 

k=m 

where 

'k' 

m 



Sl{n,k)q-\ 



[k\[k-l\---[k-m + l] 
[to] [to — 1] • • • [1] 
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is the (/-binomial coefficients. Our approach gives the recurrence relation 

L(n, m) + ^ - q J -L(n, m + 1) - qL(n + 1, m + 1) = 0. 



Similarly, for the sum 



L(n,m) = ^(-lf- fc ^,m)(T fc , 

k=0 L Jg 

we have 



(1 - g n+1 )L(n, m) + L(n + 1, m) + ^ x_ g L ( n ' m + ^ 



Exponential functions 

Noting that 



/ m+l _ 2 , a _ 1) 

+ g _ i L ( n + 1,^ + 1) — qL(n + 2, m + 1) = 0. 



A;™ = res ■ ' 



x (1 — A;x)x n+1 ' 

we can use the residue method to deal with sums involving k n . For example, consider the sum 

Applying the extended Zeilberger's algorithm to the summand 

' m \ ( — l)m-k 

F(n,m,k) - 1 

we find that 



kj (1 - fcr)x n+1 ' 



— (m + l)L(n, m) — (m + l)L(n, m + 1) + £(h + 1, m + 1) = 0. 

Since m!S2(n, m) satisfies the same recurrence relation and has the same initial values, we 
finally derive that (see [9, identity (6.19)]) 



£^fc"(-ir- fc = m!S 2 (n, 



m) 



Bernoulli polynomials 



Identities involving Bernoulli and Euler numbers have been verified in [4]. Here we only 
point out that we may also use the extended Zeilberger's algorithm to derive differential equa- 
tions satisfied by the sum. We take the Bernoulli polynomial B n (x) as an example. Recall the 
generating function 



OO t y 

n=0 



We have 



L(n, X ,y) = £ (l)y n - k B k (x) = res £ Qy" 
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The extended Zeilberger's algorithm generates 

9 , s , . j . ( (n + l)(n + 1 — k — ty) 

— F(n + 1, k, x, y) - (n + l)F(n, fc, x, y) = A fe ^ + 1 - fc) F(n ' fc ' X ' y) 

We thus have 

d 

— L(n + 1, x, y) = (n+ l)L(n, x, y). 
This relation together with the fact L(n,0,y) = B n (y) indicates that 

Y,(^)y n ~ k B k {x) = B n {x + y). 



6 Applicability of the residue method 

We have shown in Section 3 that for sums involving Stirling numbers, the residue method 
succeeds if the creative telescoping algorithm works whereas the converse is uncertain. In this 
section, we consider sums of the form 

^2F(n,k)a k , 

k 

where F(n, k) is a hypergeometric term and the generating function of at is independent of 
k. By the residue method, we aim to find a finite set S and (k, z)-free polynomial coefficients 
{Pa(ra)}cteS sucn that 

2^ P<*y n ) jfc+I = A kG(n, k, z). 

We will show that in most cases, the above equation holds only for G(n, k, z) = 0. In this case, 
we have 

^2 Pa (n)F(n + a,k) = 0, 

aeS 

which is exactly the equation appears in Sister-Celine's method. 

We first give a lemma on the C-finiteness of hypergeometric terms. 

Lemma 6. Let f(k) be a hypergeometric term and 

/(fc+i) _ m c(k+i) 

f(k) B(k) C(k) { ' ; 

be the GP -representation ((see [16] for the definition). If f(k) is C -finite, then 

A(k) = B(k) = 1. 

Proof. Suppose that f(k) is C-finite, this is, there exist constants ao, a\, . . . , ad, not all zeros, 
such that 

aof(k) + ai f{k + 1) + • • • + a d f(k + d) = 0. 
Dividing f(k) on both sides and substituting (6.1), we derive that 

A(k)C(k + l) 2 A(k)A(k + 1) C(k + 1) 
00 + aiUZ B(kj~cW~ + ^ Z B(k)B(k + l) C(k) + • • • 

.„ vd UtoMk + i) C(k + i) _ 



lfifl(H0 C(k) 
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Hence, 

d i-l d-l 

j2 ^u i c{k + i) n A(k + j) n + j) = o. 

i=0 j=0 j=i 

Since A(k) divides all the terms of the left hand side except the first one, it must also divides 
the first term. By the definition of GP-representation, A{k) is co-prime to C{k) and B(k + j). 
We thus deduce that A{k) = 1. With a similar discussion, we derive that B(k) = 1. I 

Now we are ready to give the main theorem. 

Theorem 7. Let n = (m, . . . , n r ) and F(n, k) be a hypergeometric term. Suppose that there 
are a finite set S QW and (k,z)-free polynomial coefficients {p a (n)}a<=s such that 

E, , F(n + a, k) . , , ,„ 

P«(n) V fe+1 ' ^ = A fc G(n, fc, z). (6.2) 

aes z 

Let 

9(k) = ^2 p a (n)F(n + a,k), 

aeS 

and 

g(fc + l) A(k)C(k + l) 
~gJkr~ U B(k)C(k) 

be the GP-representation. Then we have A(k) = B(k) = 1. 

Proof. Since F(n, k) is hypergeometric, there exists a rational function R(k, z) (since n is 
irrelevant, we omit these variables) such that 

G(n, k, z) = R(k, z)F(n, k, z). 

Multiplying both sides of (6.2) by z k+l /F(n,k), we see that 

h(k) = r ( fc) fi ( fc + M) _ R ^ z) 

is independent of z, where 

rW ._F(n,t + l) 



F(n, k) ' 

Suppose that R(k,z) = P(k, z)/Q(k, z), where P(k,z) and Q(k,z) are relatively prime poly- 
nomials in k, z. Then 

r(k)P(k + 1, z)Q(k, z) - zP(k, z)Q(k + l,z) = h(k)zQ(k + 1, z)Q(k, z). (6.3) 

Noting that r(k) and h{k) are independent of z, by comparing the degrees in z of both sides, 
we obtain deg^ P(k, z) = deg 2 Q(k, z). 

We first prove that z \ Q(k,z). Suppose on the contrary that there is a positive integer m 
such that z m | Q(k,z) but z m+1 \ Q(k,z). By (6.3), we see that z m+1 \ P{k + l,z)Q(k,z). 
Therefore, z \ P(k + 1, z) and hence, z \ P(k, z). But this contradicts the condition that P(k, z) 
and Q(k,z) are relatively prime. 

Then we show that Q(k, z) is independent of k. For any irreducible factor p(k, z) of Q(k, z), 
we deduce from (6.3) that p(k, z) \ zP(k, z)Q(k+l, z). Since z \ Q(k, z) and P(k, z),Q(k, z) are 
relatively prime, we have p(k, z) \ Q{k + 1, z), which implies p(k — 1, z) \ Q(k, z). By iterating 
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the above discussion, we get p(k — i,z) | Q(k, z) for any nonnegative integer i. Therefore p(k, z) 
must be independent of k. Since p(k, z) is an arbitrary factor of Q(k, z), we obtain that Q(k, z) 
is independent of k. 

Prom (6.3), we see that z \ P(k + 1, z). So we assume that 

d d+1 

P(k, z) = z^ Pi {k)z\ Q(k, z) = Y^ ft**. 

i=0 i=0 

where all qi are independent of k. Substituting these expressions into (6.3) and comparing the 
coefficient of each power of z, we find that 

r(k)p (k + l) = q h(k), (6.4) 
r(k)p 1 {k + l)-p {k)=q 1 h(k), (6.5) 

: (6.6) 

-p d {k) = q d+1 h{k). (6.7) 
By (6.4), we have po(k + 1) = qoh(k)/r(k). Substituting it into (6.5), we get 

h(k + l) h(k + l) 

Pl {k + 2) = q \ / + qi - V 



r(fc)r(fc + l) ^ r(k) ' 
Continuing this discussion, we finally derive that 

q h(k) + qi r{k)h{k+l) + q 2 r{k)r{k+l)h(k+2) + - ■ ■ + q d+1 r(k)r(k+l) ■ ■ ■ r(k+d)h(k+d+l) = 0. 
Thus the hyper geometric term 

fc-i 

v(k) = h(k)Y[r(i) 

is C-finite. Clearly, 

v(k + l) _ h(k + l) _ g(k + l) 



v(k) h{k) v ' g{k) 

By Lemma 6., we deduce that A(k) = B(k) = 1. 
For example, we consider the sum 



L(n,m) = ^2 [l) Bm+k ' 

h—n ^ ' 



k=0 

where B n is the Bernoulli number defined by 



EBn—r - — r- 
n! e z — 1 

n=0 



Rewrite the sum as 

n+m 



L(n,m)=Y J [ k " m )B k 



k=m 

Denote the inner summand by F(n, m, k). Take 



S = {(i,j):0<i,j < 1}, 
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and denote 

9(k) = ^2 Pi,j{n,m)F(n + i,m + j,k). 
(i,j)es 

We find that 

g(k + 1) _ k-m — n — 2 P(k + 1) 
g{k) ~ k + l-m P{k) ' 

where -P(fc) is a certain polynomial in k. By Theorem (6.2), we must have g(k) = 0. There is 
a non-trivial solution 

Po,o = Po,i = ~Pi,o, 

which implies that 

L(n, m) + L(n, m + 1) — L(n + 1, m) = 0. 

This coincides with the recurrence relation given by Chen and Sun [4], wherein all identities 
involving only one Bernoulli number are of this case. 
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